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ABSTRACT 

We present evidence that the 8 /jm (dust) and 24 /im luminosities of star-forming galaxies are both 
strongly correlated with their 1.4 GHz and Ha luminosities over a range in luminosity of two-to-three 
orders of magnitude. At the bright end, the correlations are found to be essentially linear over a 
luminosity range of about two orders of magnitude (corresponding to star-formation rates of several- 
tenths to several tens of solar masses per year). However, at the faint end, there appears to be a 
slope change for dwarf galaxies, possibly due to the lower dust-to-gas ratios and lower metallicities 
of the dwarfs. The correlations suggest that PAH features and mid-IR continuum emissions are good 
measures of the star formation rates of galaxies, and we present calibrations of star-formation rates 
based on existing radio and Ha relations. Our findings are based on a sample of star-forming galaxies 
selected from the main field of the Spitzer First Look Survey with the aid of spectroscopic data from 
the Sloan Digital Sky Survey and VLA 1.4 GHz data. 

Subject headings: galaxies: starburst - infrared: galaxies - stars: formation 



1. INTRODUCTION 

Star formation rates (SFRs) are fundamental measures 
of galaxy formation and evolution. Reliable measure- 
ment of SFRs could provide us with important clues and 
constraints concerning the evolution of galaxies. With 
the less extinction and the less contamination from cir- 
rus, the mid-infrared (MIR) might be ideal tracer of star 
formation. In addition, broad emission MIR features, 
such as 3.3, 6.2, 7.7, 8.6, 11.2, 12.7, 16.4, 17.1 /im, have 
already been found to be present in most star-forming 
galaxies by both the Infrared Space Observatory (ISO, 
Kessler et al. 1996; Sturm et al. 2000) and Spitzer Space 
Telescope (Smith et al. 2004), and these features are 
thought to be related to the vibrational emissions of poly- 
cyclic aromatic hydrocarbons (PAHs). 

The ISO 6.75 /im band just covers the PAH features 
whilst the 15 /im band just covers the dust continuum 
from very small grains (VSGs). However, it turns out 
that PAH emission from the 12.7 /im feature can also 
appear in the 15 /im band, so observed 15 /jm fluxes will 
not be due to VSG continua exclusively. Elbaz et al. 
(2002) have demonstrated a relationship between MIR 
(both 6.75 /im and 15 /im) and IR luminosities; whilst 
Roussel et el. (2001) and Forster-Schreiber et al. (2004) 
have demonstrated one between MIR (both 6.75 /im and 
15 /im) and Ha luminosities in spiral disks. 

With the launch of the Spitzer Space Telescope 
(Werner et al. 2004), it become available to determine 
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SFRs of millions of galaxies in the MIR bands. The 
IRAC 8/tm band covers the strongest 7.7 /im PAH fea- 
ture, whilst the MIPS 24/im band avoids most of the 
PAH features and can sample the continuum emission of 
VSGs better than the ISO 15/im band can. Calzetti et al. 
(2005) first obtained the tight relations between Spitzer 
MIR (8 /im and 24 /im) and Pa luminosities based on the 
new Spitzer observations of the 42 knots on M51. The 
goal of our work is to examine whether Spitzer Space 
Telescope IRAC 8 /im and MIPS 24/im luminosities can 
be used to measure SFRs of whole galaxies through com- 
parisons with known SFR indicators such as Ha and ra- 
dio luminosities. However, since the MIR continuum can 
easily be contaminated by emissions from dust heated 
by AGNs, special care is needed to avoid galaxies with 
AGNs. 

In Section 2, we describe how we define our sample of 
star-forming galaxies. Correlation analysis between MIR 
luminosities and radio/Ha luminosities, and calibrations 
of MIR luminosities to SFRs are presented in Section 3 
and 4. Throughout this paper, we adopt a cosmology 
with n M = 0.3, Q A = 0.7, and H = 70. 

2. SAMPLE SELECTION 

Our galaxy sample was taken from the extragalactic 
component of the main field of the Spitzer First Look 
Survey (FLS). In particular, it was drawn from the over- 
lap area of about 3.7 degree 2 of the main field that has 
been imaged by both instruments: IRAC (3.6, 4.5, 5.8, 
8.0 /im) and MIPS (24/im). This area has also been cov- 
ered by the SDSS (Stoughton et al. 2002). The IRAC 
images (in all four IRAC bands) were mosaiced from the 
"Basic Calibrated Data" (BCD) as described by Fazio, 
et al. (2004) and Huang et al. (2004); whilst the MIPS 
images were mosaiced from post-BCD images. Matching 
the sources detected by SExtractor (Bertin & Arnouts 
1996) in all five bands with the 2mass sources ensured 
astrometric uncertainties of less than 0.1". The IRAC 
and MIPS bands have flux calibration accuracies of bet- 
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ter than 10% (Fazio et al. 2004; Rieke et al. 2004). 

FLS sources were then matched with the main galaxy 
sample of the second data release of the SDSS (Strauss et 
al. 2002) with a matching radius of 2". The galaxies with 
Ha fluxes of S/N > 5 were selected. Based on the tradi- 
tional line-diagnostic diagram [Nil] /Ha versus [OIIIJ/H/3 
(Veilleux and Osterbrock 1987; Wu et al. 1998), narrow- 
line AGNs were removed. The fluxes of all of the emission 
lines were taken from the SDSS catalogue of measured 
line fluxes (Tremonti et al. 2004, version 5.0.4). We also 
included 1.4 GHz fluxes when the coordinates of VLA 
1.4 GHz radio catalogue of Condon et al. (2003) matched 
the sample galaxies to an accuracy of 0.5 ". 

Our sample of the star-forming galaxies contains a to- 
tal of 91 galaxies with both Ha and MIR fluxes in all four 
IRAC bands. 70 of these also have MIPS 24/im fluxes. 
Due to sensitivity limitations affecting the VLA obser- 
vations (Condon et al. 2003), only 38 of our 91 IRAC 
sources and 33 of our 70 MIPS sources have 1.4 GHz 
fluxes. All these objects are local with redshifts of less 
than 0.2. There are 12 dwarf galaxies with absolute B 
magnitudes fainter than —18.0 in our sample, but only 
three of these dwarfs have radio data. 

3. MIR LUMINOSITIES AND CORRELATION ANALYSIS 

To examine whether MIR 8/im and 24/mi luminosities 
can be taken as measures of the SFRs of star-forming 
galaxies, we first compare them with 1.4 GHz and Ha 
luminosities, which have already been extensively used 
as measures of SFRs. To correct for redshift effects, we 
adopt a steep non-thermal radio spectrum whose flux 
obeys a power law of the form v~ a (Cardiel et al. 2003) 
and has a spectral index of a = 0.8 (Condon 1992), and 
compute rest-frame radio luminosities at 1.4 GHz. 

As for the Ha luminosities, both extinction and aper- 
ture corrections were applied to derive the total Ha emis- 
sion of the whole galaxy based on the fibre spectra. The 
extinction correction included both the foreground ex- 
tinction of the Milky Way and the intrinsic extinction of 
galaxy itself. The Milky Way extinction was corrected by 
the parameterized curve of Cardelli et al. (1989) and the 
intrinsic extinction An a was obtained from the Balmer 
decrement Fho/Fh/s (Calzetti 2001). For the aperture 
corrections, we followed the Formula Al of Hopkins et 
al.'s (2003). The corrections were based on the r band 
fiber magnitude and the corresponding Petrosian mag- 
nitude, which can be approximately treated as the total 
flux of the whole galaxies in r band. 

Using SExtractor, "AUTO" magnitudes were obtained 
for the total luminosities of the sample galaxies in the 
IRAC and MIPS 24/im bands. We adopted the SED 
of the normal HII galaxy NGC 3351 from IRS observa- 
tion of the Spitzer Legacy Program SINGS (Kennicutt 
et al. 2003) as the template for the K-corrections for our 
sample galaxies in both the 24/im and 8/im bands. To 
estimate the uncertainty of SED template, we compared 
it with those of other 20 SINGS normal galaxies, M82, 
and the models of normal and starburst galaxies of La- 
gache, Dole and Puget (2003). Within redshift of 0.2, 
the difference of K-correction is at most 30%. Although 
PAH emissions dominate the 8.0/im band for most of the 
sample galaxies, there is still a stellar continuum in this 
band, especially for the galaxies of earlier type. To re- 
move the stellar contribution in this band, we made use 



of the 3.6/im band which avoids the strong PAH emis- 
sion and used it to estimate the stellar contribution in 
the 8.0/im band. We selected 74 galaxies without optical 
emission lines in the FLS and assume they are dust-free 
at 3.6/im. The statistics show that the averaged mag- 
nitude differences between the 3.6/im and 8.0/im bands 
for these galaxies is about 1.45 mag with a scatter of 
0.28 mag for z < 0.2. This corresponds the factor of 
0.26 to scale the stellar continuum of 3.6/im to that of 
8/im, and a little larger than the 0.232 of Helou et al. 
(2004). This factor was then used to subtract the stellar 
continuum contributions from the estimated luminosities 
for sample galaxies in the 8/im band. Generally, the 
stellar continuum flux amounted to about 10% of the 
8.0/im flux, similar to that of Engelbracht et al. (2005). 
The resulting luminosities which represent the dust emis- 
sions of our sample galaxies will henceforth be denoted 
8/im(dust). Though the PAH emissions dominate the 
8/xm band, but the dust grains with size > 50A can still 
contribute to the 8/im continuum at some level (Li & 
Draine 2001). 

Figure shows (a) 24/im and (b) 8/im(dust) luminos- 
ity versus the 1.4 GHz radio luminosity in the log-log 
space. We can see that a good correlation exists between 
the MIR and radio luminosities over a luminosity range 
of about 2 orders of magnitudes. The unfilled symbols 
represent the only three dwarf galaxies with catalogued 
radio fluxes. These three dwarfs seem to follow similar 
correlations to those found for brighter galaxies. The 
solid lines in the figures are the best fits to all brighter 
galaxies with two-variable regression. The dotted lines 
are the linear fits. The fitted parameters and the corre- 
lation coefficients are listed in Table 1. As evident from 
Figure GJa), our results are consistent with those of Ap- 
pleton et al. (2004). However, our correlation exhibits a 
greatly reduced scatter of 0.17 cf. Appleton et al.'s (2004) 
0.27. In the case of Figure^(b), our results appear to be 
consistent with those of Elbaz et al. (2002) at the faint 
end but there is a departure at the bright end. The corre- 
lation was deduced from the Elbaz et al.'s relation of the 
bolometric IR luminosity versus ISO 6.75/im luminosity 
and the radio "q" parameter. 

FigureElshows (a) 24/im and (b) 8/im(dust) luminosity 
versus Ha luminosity in log-log space. We can observe 
a good correlation between these two MIR luminosities 
and the Ha luminosity for the sample galaxies except 
that there appears to be a slope change for dwarf galax- 
ies. Both of the linear correlations for non-dwarf galax- 
ies have correlation coefficients of about 0.9. As evident 
from Figure (b) broad agreement is found with the 
results of Roussel et al. (2001) which were based on a 
straight-line fit to 6.75/im (ISO) and Ha luminosities for 
spiral disks. 

4. CALIBRATION OF 24/iM AND 8//M(DUST) STAR 
FORMATION RATES 

From the strengths of the correlations listed in Table 1, 
there appears to be a strong link between both MIR 
24/im and 8/im(dust) luminosity and both 1.4 GHz radio 
and Ha luminosity. As both radio and Ha luminosities 
are already considered to be good measures of SFR, we 
can estimate SFRs based on MIR (24/im and 8/im(dust)) 
luminosities from either radio and/or Ha SFRs. 

Combined with the linear relations in Table 1 and the 
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SFR-radio luminosity relation given by Yun, Reddy and 
Condon (2001), we can derive SFRs from both the 24/im 
and 8/im(dust) luminosities: 



SFR 2 ^ m {MQyr- 1 ) = 



SFR Slxm ( duat )(M Q yr x ) 



v L u [2A/i to] 
6.66 x 1O 8 L ' 

v L u [8fj,m(dust)] 
1.39 x 10 9 L m ' 



(1) 



(2) 



Similarly, we can also derive SFRs from Ha luminosi- 
ties based on Kennicutt's (1998) SFR formula: 



SFR 24tim (M Q y r - 1 ) = 



SFR^ m ( dust) (M Q yr r ) 



v L„[24// to] 
6.43 x 1O 8 L ' 

v L l/ [8/im(dust)] 
1.57 x 10 9 L m ' 



(3) 



(4) 



However, due to the change in slope exhibited by the 
dwarf galaxies in Figure |3 the latter two relations are 
necessarily limited to the normal star-forming galaxies. 

There are many factors which could affect calibrations 
and increase the scatter in the correlations. These in- 
clude the accuracy of the fibre-aperture corrections, the 
validity of using the Balmer decrement to estimate ob- 
scuration in galaxies (Hopkins et al. 2003) , possible radio 
and MIR contamination from obscured weak non-optical 
AGNs, the template used for K-correction, and so on. We 
have therefore compared the SFRs derived from formulae 
(1) through (4). Reassuringly, we found that the agree- 
ment between the SFRs based on radio and Ha calibra- 
tions are around 10% for both the 24/xm and 8/im(dust) 
relations respectively. This level of error is well within 
the 1-er scatters as listed in table 1. 

5. DISCUSSION 

Based on the relations we have derived, we suggest 
that the 24/im continuum could be due to the warm 
dust component associated with star formation regions 
(O'Halloran et al. 2005). This support the case made by 
Dale et al. (2001) that the 20-42 /im continuum may be 
the best dust emission tracer of current star formation in 
galaxies. The results of the HII knots in M51 (Calzetti et 
al. 2005) and the 75 nearby galaxies from SINGS (Dale 
et al. 2005) also suggest that the 24/im emission is a 
useful tracer of the local star formation. Furthermore, 
the agreement between our radio-24/im relation for star- 
forming galaxies and Appleton et al.'s (2004) relation for 
radio-quiet sources in the same Spitzer FLS region (Fig- 
ure H (a)), indicates that even for many optical AGNs, 
the 24/im emission is still dominated by star formation. 
The much larger scatter in the radio-24/im relation for 
radio-quiet sources could then be explained by the extra 
contribution from the central AGN. 

Though we have drawn the linear calibration from MIR 
luminosity to SFRs, the best fit slopes for the 8/im(dust) 
luminosity versus Ha/radio luminosity is obviously shal- 
lower than those for 24/im (Table 1). Calzetti et al. 



(2005) obtained a similar results for HII knots in M51. 
Their 24/im and Pa luminosities are nearly linear corre- 
lated, but the slope for 8/xm luminosity versus Pa lumi- 
nosity is much shallower than one. Peeters, Spoon and 
Tielens (2004) have pointed out that PAHs may be bet- 
ter tracers of B stars than tracers of massive star forma- 
tion (O stars). PAH features can be strong in moderate 
UV radiation fields such as photo-dissociation regions, 
but the PAH carriers could be destroyed in hard radia- 
tion field (Galliano et al. 2005). If it is true, 8/im(dust) 
flux might not be a linear measure of intensive star- 
burst activity, this could explain the shallower slope of 
8/im(dust). Alternatively, the VSG continuum might be 
a useful tracer of star formation in higher UV-radiation 
fields (Forster-Schreiber et al. 2004), though VSGs might 
also be destroyed in very harsh radiation fields (Contursi 
et al. 2000). Since the sample used here are star-forming 
galaxies, a sample of luminous or even ultraluminous in- 
frared galaxies is needed to check if such relations could 
be extend to intensive starburst environments. 

Kewley et al. (2002) showed that SFRs derived from 
Ha emission are similar to those derived from IR emission 
even at SFR of 0.01 M yr _1 . However Figure |21 shows 
that the dwarf galaxies systematically deviate from the 
correlation of the star-forming galaxies to the side of 
weak PAH emission in the log-log space. Madden (2000) 
and Dale et al. (2005) have also shown that PAH fea- 
tures are weak in dwarf galaxies. It is well known that 
many dwarf galaxies are gas-rich and have higher gas-to- 
dust ratios (as well as lower metallicities) than normal 
galaxies. In low metallicity environments, the UV radi- 
ation field is presumably strong enough to destroy the 
PAH carriers and weaken the PAH features, for exam- 
ple, Blue compact dwarf galaxies SBS 0335-052 (Houck 
et al. 2004). Estimating the oxygen-abundance with the 
emission line ratios (Kobulnicky, Kcnnicutt, & Pizagno 
1999), we find that two of four dwarf galaxies (including 
the most deviant one) exhibit much lower metallicities 
(about one half or one-third of the solar metallicity) than 
those of the normal star-forming galaxies. So the slope 
of the MIR-SFR relations might therefore be metallicity 
dependent. This is also supported by Engelbracht et al. 
(2005) 's result that the 8/im-to-24/im flux ratio strongly 
depends on metallicity for a sample of 34 star-forming 
galaxies. 
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20 21 22 23 

Logio Li.4ghz (W Hz" 1 ) 

Fig. 1. — The correlations between the MIR and radio luminosities. The filled circles are the sample galaxies and the unfilled symbols are 
the dwarf galaxies with Mg < —18.0 mag. The solid and dotted lines are the best non-linear and linear fits respectively for star-forming 
galaxies, (a) 24/im versus radio (1.4 GHz). The long-dashed line represents the linear fit to the radio and 24/im luminosities obtained by 
Appleton et al. (2004) for radio-quiet sources with K-corrections based on an M82 template, (b) 8/^m(dust) versus radio. The short-dashed 
line represents the relation of the radio and 6.75/im (ISO) luminosities drawn from Elbaz et al. (2002). 
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6 7 8 9 

Log 10 L Ha /Lo 

Fig. 2. — The correlations between the MIR and Ha luminosities. The filled, unfilled circles, the solid and dotted lines are as in Figure 1. 
(a) 24^tm versus Ha. (b) 8/im(dust) versus Ha. The dot-dashed line is the linear relation between ISO 6.75/im and Ha luminosities for 
spiral disks found by Roussel et al.s (2001). 
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Table 1 

Correlation coefficients for luminosity relations for star-forming galaxies 



y 


X 


a 


b 


r 


c 


N 


vL v \p.Ay,m](Lo) 


L 1A gh z (W Hz~ L ) 


-13.31±0.42 


1.04±0.09 


0.82 


-12.41±0.17 


30 


v L„[8/im](Z/0) 


L 1A GHz(W Hz' 1 ) 


-10.12±0.31 


0.91±0.07 


0.88 


-12.09±0.13 


35 


vL v [24fxm](LQ) 


Lhc{Lq) 


0.34±0.19 


1.12±0.07 


0.88 


1.29±0.19 


63 


uL v [8fxm](LQ) 


Lh<x{Lq) 


2.31±0.15 


0.92±0.05 


0.87 


1.68±0.19 


79 



a and b are the coefficients of the non-linear fitting of log 10 (y) = a + blog 10 (x) and c is the coefficient of the linear fitting of log 10 (y) = 
c + log 10 (a;). r is the Spearman's correlation coefficient, rho, and N is the number of sample galaxies (excluding the dwarf galaxies) used 
for the fitting. 



